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Abstract 

Photon splitting cascades in a magnetic field are considered. It is shown, in 
the stochastic limit approximation, that photons in cascades might form entangled 
states ("triphons") and that they obey not Bose but a new type of statistics, the 
so called infinite or quantum Boltzmann statistics. These states belong to an inter- 
acting Fock space which is a generalization of the ordinary Fock space. The new 
photon statistics in principle can be detected in future astrophysical experiments 
such as the planned Integral mission and also in nonlinear quantum optics. 
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One of the most interesting manifestations of the non-linearity of Maxwell's equations 
with radiative corrections is the splitting of a photon into two in an external magnetic 
field. In a constant uniform field, this process occurs with conservation of energy and 
momentum. The process was considered by Adler et al. in the early '70s by using the 
Heisenberg-Euler effective Lagrangian |T|, 0, [|. Photon splitting was considered as a 
possible mechanism for the production of linearly polarized gamma-rays in a pulsar field. 
Recently the splitting of photons has found astrophysical applications in the study of 
annihilation line suppression in gamma ray pulsars and spectral formation of gamma ray 
bursts from neutron stars @, [5J. Photon splitting cascades have also been used in models 
of soft gamma-ray repeaters, where they soften the photon spectrum || f|. The process 
of photon splitting is potentially important in applications to a possible explanation of 
the origin of high energy cosmic rays from Active Galactic Nuclei ||. A recalculation of 
the amplitude for photon splitting in a strong magnetic field has been performed recently 

m uniini. 

In this note we argue that photon splitting cascades in the magnetic field might create 
entangled states and that photons in cascades obey not Bose but a new type of statistics, 
the so called infinite or quantum Boltzmann statistics. These states are formed from 
triples of entangled photons and may be called triphons. They belong to a generalization 
of the ordinary Fock space which is called an interacting Fock space. Creation and an- 
nihilation operators for infinite (quantum Boltzmann) statistics satisfy neither Bose nor 
Fermi commutational relations but 



6(k)6 + (p)=5(k-p) (1) 

The relations (p]) have been considered in several recent works in quantum field theory. 
The second quantized example of infinite statistics has been discussed by Greenberg [|1^ 
and in the present note it is shown that this statistics has a physical meaning since it 
describes photons in cascades and more generally the dominating diagrams in the long 
time/week coupling limit in quantum field thory. The notion of interacting Fock space 



was introduced by Accardi and Lu |T3[ in nonrelativistic QED and it was related with 
the role of non-crossing diagrams in the stochastic limit. The master field describing the 
large N limit in QCD was obtained in [ 14| , it is quantized by using the relations ([Xj) ; see 



15, [L6| for a recent discussion of the large N limit. 

We will start from a discussion of the theory of photon splitting cascades and show 
the emergence of infinite statistics in this theory and then discuss its connection with the 
stochastic limit of quantum field theory. 

In the decay of a photon with momentum k into photons with momentum ki and k2, 
we have conservation of momentum and energy k = ki + k2 , w(k) = ui(k\) + 002(^.2) ■ 
For photons in vacuum, in the absence of external fields, u = u>i = u 2 = k and although 
these two equations have a solution the decay is forbidden by the invariance under charge 
conjugation (Furry's theorem). 

In a constant uniform magnetic field Bo there are only two decay processes kinemati- 
cally allowed, 7y — > 7_i_ + 7± and 711 — > 7y+7± @. Here the subscripts _L and || will denote 
polarizations of the photon with respect to the vector Bo- More precisely, in presence of 
a magnetic field one has a distictive plane, namely the kB plane. One takes the linear 
polarization of the magnetic field of the photon parallel and orthogonal to this plane as 
the two independent polarizations of the photon, || and _L, respectively. 
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The vacuum in the presence of the field Bo acquires an index of refraction n, and the 
photon dispersion relation is modified from k/uo = 1 to k/uo = n. The indices of refraction 
nil j_ can be calculated from the Heisenberg-Euler effective lagrangian. Adler showed 
that for subcritical fields in the limit of weak vacuum dispersion only the splitting mode 
|| — >_L + _L operates below pair production threshold. For weak dispersion n± = 1 + ^ f3 

and n\\ = 1 + Jr f3, where f3 = B 2 sin 2 9 and 9 is the angle between k and B . It is 
mentioned by Harding et al. that in magnetar models of soft gamma repeaters ||, where 
supercritical fields are employed, moderate vacuum dispersion arises. In such a regime, 
it is not clear whether Adler's selection rules still endure since in his analysis higher 
order contributions to the vacuum polarization are omitted. In [[J photon cascades are 
considered for the case where all three photon splittings modes allowed by CP invariance 



are operating. Baier et al. [JTOJ have found that there is only one allowed transition 
(|| — >J_ + _L) for any magnetic field. They suggested that a photon cascade could develop 
only if magnetic field changes its direction. It seems that the question on the validity 
of Adler's rule for a non weak vacuum dispersion deserves a further study. In this work 
we consider photon cascades when both kinematically allowed modes (|| + _L and 
|| — > ||+ -L) operate. 

The interaction operator for the decay || + _L is known to be || 

V 1 (t) = X 1 J (B E 1 )(B E 2 )(B B)rf 3 x, (2) 

where the coupling constant Ai = 13e 6 /3157r 2 ra 8 and magnetic and electric parts of photon 
field are 

B = z(47r) 1/2 k x eiie-^-^antki), Ex = -z(47r) 1/2 cj 1 e ± e i(kir -" 1 ' ) a+(k) (3) 
and similarly for E 2 . 

For the decay || — > 1 1 — I — I one has a similar interaction operator with the operator 

structure 

A+(t) = Aa+(k 1 )a+(k 2 )a||(k)e-^(5(k - k x - k 2 ) (4) 

where E = cjii(k) — a;ii(ki) — a;j_(k 2 ). The coupling constant A in this case can be estimated 
as A/Ai = a(B /B cr ) 2 , where a is the fine structure constant, a = e 2 /Tic and B cr = 
m 2 c 3 /eh ~ 4.4 x 10 13 Gauss. 

Let us consider a photon cascade created by a photon with momentum k and polar- 
ization || . The photon splits as 711 (k) — > 7_i_(ki) + 711 (k — kj.). Then one has the next 
generation of splitting: 7y (k — k x ) — > 7_i_(k 2 ) + j\\ (k — k x — k 2 ) etc. After generations of 
splitting one gets a cascade with A^ photons with _L polarization and momenta k 1; . . ., k^v 
and also one photon with || polarization and momentum k — ki — . . . — kyv. An example 
of a cascade with two generations is shown in Fig 1. 

It is important to notice that we consider cascades with real photons and therefore the 
diagram in Fig.l is not a Feynman one because all the lines (including an intermediate one) 
correspond to real particles on the mass shell and not to virtual states. From the point 
of view of the standard quantum field theory all the lines in the diagram are "dressed" 
lines on the mass shell and moreover the initial photon 711 (k) is prepared in a special way 
such that it undergoes the decay in a finite time. So we cannot use the standard S"-matrix 
approach and the standard Feynman diagram technique to describe this process. The 
diagram is also not a diagram in the non-covariant diagram technique [0 because we 
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Figure 1: Cascade with two generations 
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Figure 2: Crossing and non-crossing diagrams 



have the conservation of energy at every vertex. The cascade in Fig.l may be intuitively 
described by the following state 



where momentum conservation is build in creation and annigilation operators and energy 
conservation is accounted for by the factor /(p, k) = f(u)\\ (p) — uj±(k) — uj\\ (p — k)) where 
f(co) is a function with support at u = 0. As we shall see below this is not the 5- 
function but roughly speaking its square root. Indeed, the transition amplitude between 
two cascade states is given by scalar product 



(^(k'XXMO^MM)) = |/(k,k 1 ))| 2 |/(k-k 1 ,k 2 )| 2 5(k-k')5(k 1 -k' 1 ) ( 5(k 2 -k' 2 ) (6) 



Notice that in the scalar product @ only the non-crossing diagram, (Fig. 2a) contributes. 
In fact the contribution from the crossing diagram in Fig. 2b vanishes because of conser- 
vation of energy and momentum. This is a crucial point where the difference between 
our diagrams describing real particles in intermediate states and the Feynman diagrams 
having virtual particles in intermediate states is evidentiated. In the Feynman diagram 
technique the amplitude of emission of the two photons is represented by a sum of two 
diagrams differing by the order in which the two photons are emitted. Here we have only 
one diagram, Fig.l. 

Now let us observe that if in (^) we replace operators a^(kx) and a^(k 2 ) by the quan- 
tum Boltzmann operators b±(ki) and 6j_(k 2 ) satisfying the relations (1), i.e. 6j_(k)6^(p) = 
5(k — p), then it will be automatically guaranteed that only the non-crossing diagrams 
survive. Therefore it is natural to describe cascade wave functions in terms of these opera- 
tors. It is well known that standard free photons are bosons. Therefore to see the quantum 
Boltzmann statistics we have to prepare a special state depending on interaction. A nat- 
ural method, leading to this result, is suggested by the stochastic limit technique. In fact 
it is natural to expect that the cascades with physical intermediate states occur at a time 



k,k 1; k 2 )> 



a+(k - ki - k 2 )/(k - k ls k 2 )a+(k 2 )/(k, k 1 )a+(k 1 )|0) 



(5) 
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scale slower than the one occurring in the standard S- matrix approach to multiparticle 
production. Notice also that the coupling constant in the interaction term (f|) are small. 
Therefore we are precisely in the situation in which one considers long time-commulative 
effects of weak interactions. The stochastic limit captures exactly these effects in the van 
Hove limit A — > 0, t — > oo so that A 2 t ~ constant = r (new time scale) which means that 
we measure time in units of 1/A 2 where A measures the strength of the self-interaction 
(proportional to the magnitude of the magnetic field in our case). It is remarkable that 
in this limit the triplets of photons ("triphons") behave like a single new quantum field 
obeying a new statistics. 

Now let us consider the question how one can prepare a state with the new statistics for 
photons. If we would deal with the scattering at infinite time (^-matrix) we simply have 
to consider two Feynman diagrams to take into account the Bose statistics of photons. 
However in the cascade we deal with evolution in finite time and the states of photons 
7_i_(ki) and 7_i_(k 2 ) are prepared in a special way because they are emitted at time t 1 
and t 2 , respectively. Therefore, there is a reason not to add the second diagram. There 
exist a special procedure which is adequated to this situation. This is so-called stochastic 
limit technic. This limiting procedure is widely used in the consideration of the long 
time/weak coupling behaviour of quantum dynamical systems with dissipation, see for 
example p|, [19[. 

In this procedure one deals with states generated by products of rescaling interaction 
in different times in the interaction picture ^ V (-£2) in the limit of A — > 0. This is 



connected with anisotropic asymptotics, see p0| , where one deals with correlators < 
V{ti/\)...V{t n /\) > when A — > 0. The equation for the evolution operator in interaction 
picture reads 

where A is the coupling constant. In the stochastic approximation one replaces U^ x \t) to 
another operator U (t) 

U w (t)^U(X 2 t) 

where lA{t) is obtained by performing the van Hove rescaling of time t —>■ t/X 2 and taking 
the limit A — > 0. Then for the limiting evolution operator U(t) = lim^o U^(t/X 2 ) one 
gets the equation 

d J§l = -iv { tw) 

where 

V(f) = lim \v ( 4 
a^o A VA 2 

For the interaction Hamiltonian (|J) we consider the asymptotic behaviour of the col- 
lective operator A\(t) = jA(j%) and its Hermitian conjugate. We obtain 

lim At(t) = b+B + (k 1: k 2 , k)5(k - ki - k 2 ) (7) 

A — >0 

where 

B+Ckx.k^k) = 6+(k 1 )6+(k 2 )6||(k)(27r) 1 /V/2 (jE) 
E = u\\(k) - W||(ki) - u±(k 2 ). 
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and the following commutation relations take place 



b t b+ = 5(t - t') (8) 

b ± (p)b + ± (p') = 6(p-p'), b ]l (p)b{(p') = 5(p-p') (9) 
and also the following relation, which explains out notation 5 1 ^ 2 (E) 

fi(k 1 ,k 2 ,k)5+(k' 1 ,k^,k') = 27r5{E)5{k 1 -ki)5(k 2 - k' 2 )5(k - k>+(k)b,|(k) (10) 

The relation (|7|) is understood in the sense of convergence of matrix elements, 

(0\a ]l (p 1 )A*(t 1 )...A*(t n )a+(p 2 )\0) -a^o (0|6||(p 1 )S#... J B t t6jf(p 2 )|0) (11) 

where A& means A or A + 

Relations (f|) define the free or Boltzmann commutation relations. Notice the Boltz- 
mannian white noise relation (§), which makes our model particularly suitable for Monte 
Carlo simulations. The origin for arising these new commutation relations lies in the fact 



that the crossing diagrams in the computation of the matrix element ( TT]) are supressed in 
the weak coupling /large time limit. The presence of the 5(i£)-factor has two important 
physical consequences. First, the commutation relations for the £?# are not a consequence 
of the corresponding relations for bj and &_[: the three photons are entangled into a single 
new object (triphon). Second, the triphon creation and annihilation operators 5* operate 
not on the usual Fock space but in interacting Fock space . 

Let us illustrate on the example of four triphons that only diagrams with non-crossing 
lines survive in the limit A —>■ 0. An arbitrary diagram with four "triphons" schematically 
can be written as 

1 / e 3 v{^ y £(±)E®t i }<l>{t,p,k)dkI[dt i 

Here are the same as in (TJ) with fc® being momenta of line coming in and ki \ 
momenta of lines coming out from z-vertex; ± correspond to vertex (|j) and its complex 
conjugated; p are external momenta and k are two independent momenta; accomulates 
all form-factors and test functions. In general, the sets of momenta corresponding to 
different verticies are different. But if there are two verticies such that momenta coming 
in the first vertex come out from the second one and via versa we call these verticies 
as conjugated ones (see Fig. 2a where conjugated verticies are denoted by hat). Only 
diagrams consistent of pairs of conjugated vertices survive in the limit A — * 0. Indeed, 
making a change of variables (ti, £ 2 , t%i t§) —> (ti, t 2 , tj, t%), T\ = t -^r L , t 2 = -^A we see 
that diagram Fig. 2a gives contributions containing the following factors 

5(^ 1 l ||+± )5( J E | fi ||+± ) ( 5(t 1 - ti)tf(<a - t^, 

where the energy factors E^\ E^ are not independent due to momentum conservation, 
that is typical for interacting Fock space |TB |. For all others diagrams a similar change 



of variables do not remove the dependence of exponent of A that due to fast oscillations 
produces zero contributions. 

A photon splits into two not only in a magnetic field but also in a nonlinear medium. 
In fact such processes are well known in nonlinear quantum optics, see for example pi] .In 



the nonlinear process of parametric down conversion a high frequency photon splits into 
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two photons with frequencies such that their sum equals that of the high-energy photon. 
The two photons produced in this process possess quantum correlations and have identical 
intensity fluctuations. The photon pairs generated in parametric down conversion carry 
quantum correlations of the Einstein-Podolsky-Rosen type. Experiments to test Bell 
inequalities were designed using a correlated pair of photons. A two photon cascade 
was used in the initial experiment by A. Aspect to generate the correlated photons but 
more recent experiments have used parametric dawn conversion. These experiments have 
consistently given predictions of quantum theory and in violation of "realistic" classical 
predictions. In our case the quantum correlations do not come from superpositions of 
spins or polarizations but have a deeper dynamical origin expressed by the 5(E) -function 
in formula (0). It would be very interesting to extend these experiments to observe the 
new statistics considered in this paper. 

In conclusion, in this note we have argued that photon cascades in a strong magnetic 
field might create entangled states (triphons) which obey not Bose but the quantum 
Boltzmann statistics. This prediction is based on the assumption that both kinematically 
allowed photon splitting modes operate. Another assumption is that intermediate photons 
in a cascade are physical particles (i.e. on the mass shell) but not virtual states. This 
is equivalent to the validity of the stochastic limit approximation. In fact both of these 
assumpions deserve a further study. A better theoretical understanding of the photon 
splitting with a non weak dispersion is required. From the experimental side new more 
precise devices such as the planned Integral misson [|J might significantly advance our 
understanding of the fundamental problem of photon statistics. 
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